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We discuss a new family of multi-quanta bound states in the Standard Model, which exist due to 
the mutual Higgs-based attraction of the heaviest members of the SM, namely, gauge quanta W, Z 

and (anti)top quarks, f, t. We use a self-consistent mean-field approximation, up to a rather large 
particle number N . In this paper we do not focus on weakly-bound, non-relativistic bound states, 
but rather on "bags" in which the Higgs VEV is significantly modified/depleted. The minimal 
number A'^ above which such states appear strongly depends on the ratio of the Higgs mass to the 
masses of W,Z,t,t: For a light Higgs mass tuh ~ 50GeV bound states start from A'' ~ 0(10), 
but for a "realistic" Higgs mass, ruH ~ 100 GeV^, one finds metastable/bound W,Z bags only for 
N ~ 0(1000). We also found that in the latter case pure top bags disappear for all N, although 
top quarks can still be well bound to the W-bags. Anticipating cosmological applications (discussed 
in a companion paper) of these bags as "doorway states" for baryosynthesis , we also consider the 
existence of such metastable bags at finite temperatures, when SM parameters such as Higgs, gauge 
and top masses are significantly modified. 



I. INTRODUCTION 
A. Motivation 

In the Standard Model (SM), the interac- 
tion of particles includes an attractive Higgs 
exchange. For a two-particle system it is not 
difficult to see under which condition a Higgs 
exchange would lead to bound states of such 
particles. Unfortunately, one finds that the cor- 
responding critical Higgs mass lies far below the 
current experimental bound > 116 GeV. 

But one should not be discouraged too early 
by this example. Being a scalar, the Higgs 
generates universal attraction between all kinds 
of particles. Furthermore, the strength of the 
attraction is proportional to their mass, sim- 
ilar in this respect to the gravitational force. 
Gravity, feeble as it is, is able to hold together 
planets, stars and even create closed systems 
(black holes), because the rather weak cou- 
pling can be compensated by a large number 
N of participating particles. Unlike vector-field 
based forces induced by electric, weak or color 
charges, both gravity and scalar exchanges are 
exempt from "charge screening" and become in- 



creasingly stronger for large number of parti- 
cles. However, there is an important difference 
with gravity in that the Higgs boson is neither 
massless, nor particularly light in comparison 
to W, Z or t. This leads to the following ques- 
tion: What happens with heavy multi-quanta 
states when the Higgs mass is increased, from 
a near-zero value, to Mh ^ 0(100 GeV) where 
it may be soon found. This is the main subject 
to be addressed in this work. 

An instructive analogy is provided by nuclear 
physics. It is convenient to describe such situa- 
tion by a (much-simplified) Walecka model, in 
which the nuclear forces can be approximately 
described by the a and co meson exchanges. 
The correlated two-pion state a is "the Higgs 
boson of the nuclear physics", obtaining VEV 
in chiral symmetry breaking. The a and u) me- 
son masses set a scale l/m^.o; ^ 0.3 /m for the 
range of nuclear forces. Naively it is too small 
compared to nuclear sizes (several fm) or to the 
typical inter-nuclear distances n~^^^ ~ 1.5 fm. 
Furthermore, because of similarity of masses 
mcr ~ 600 MeV, ~ 770 MeV, as well as cou- 
plings, the sigma-induccid attraction is nearly 
canceled by the omega-induced repulsion. The 
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sum is an order of magnitude smaller than one 
would get from scalar and vector components 
taken separately. The resulting nuclear forces 
thus miss to bind two neutrons and can barely 
bind a deuteron. And yet there are lots of 
bound/metastable nuclei, a subject of active 
branch of physics for the last century. 

Is the situation at electroweak scale similar? 
It is widely believed that the Standard Model 
is just a low energy effective Lagrangian, hiding 
important physics behind its simplistic scalar 
Higgs. Indeed, both in superconductivity and 
in QCD chiral symmetry breaking, there are no 
fundamental scalars, and the condensates are 
in fact the two fermion or fermion-antifermion 
combinations, respectively. The fermionic mass 
generation by the Yukawa couplings is perhaps 
also just a parameterization, soon to reveal 
its dynamical content in coming LHC experi- 
ments. We are not going to speculate on this is- 
sue which has huge literature, mentioning only 
one thing. We are not aware of any particular 
model which would proposes a vector compan- 
ion to Higgs with a similarly small 0(100 GeV) 
mass. For example, the "techni-p" is predicted 
to be at the scale Mp = 1 — 2TeV, see e.g. 
[I]). Thus, unlike in the nuclear physics, one is 
not expecting the scalar-vector cancellation to 
take place. Heavier exchanges with a TeV scale 
mass can be included as perturbation later, and 
for now effective description of the SM Higgs 
seems to be quite adequate for the purposes of 
this paper. 

In this work we study the binding of vari- 
ous multi-quanta states within the SM, both 
semi- analytically for a large number of particles 
N ^ oo, and numerically for finite N. We will 
discuss the conditions under which metastable 
and bound states of N heavy quanta occurs, 
considering heavy fermions (top quarks) and 
gauge bosons {W, Z). We study these objects in 
vacuum (zero temperature/density), as well as 
at finite temperature, envisioning possible cos- 
mological applications. The method used is the 
mean-field approximation, ignoring 0{1/N) ef- 
fects such as center of mass motion. Unlike pre- 
vious works, we do not focus on weakly bound 
nonrelativistic systems but rather on "bags" in 
which the Higgs VEV inside is strongly modi- 
fied. (One strong motivation for that is related 
with occurrence of the electroweak sphalerons 
inside them, see the companion paper re- 



quiring near-empty bags.) We will see that bags 
containing W-bosons always exist, regardless of 
the value of the Higgs mass, provided there's a 
large enough number of them. Next, we study 
pure top bags and find that for a wide range of 
N pure (anti-)top bags at realistic Higgs mass 
are excluded, even as metastable minima. 

Why should one study these multi-quanta 
states? From a methodical point of view, they 
are a new class of manybody systems, beyond 
atoms and nuclei, which are truly relativis- 
tic. Although the experimental production of 
a large number of heavy quarks/bosons in a 
sufficiently small volume is clearly impossible, 
at the LHC or any other proton accelerator , 
their existence may be important cosmologi- 
cally. In fact, in our companion paper [3] we 
show that such bound states may play a sig- 
nificant role in cosmological Baryogenesis, pro- 
vided they significantly deplete the Higgs VEV 
and are sufficiently long-lived. We specifically 
consider a scenario, studied theoretically and 
numerically for some time, in which the transi- 
tion to the current broken phase of the SM hap- 
pens directly from inflation reheating, avoiding 
the electroweak phase transition. 



B. Recent works 

Having outlined the main issues to be ad- 
dressed in this paper, let us provide an overview 
on recent related studies. The existence of bags 
due to a Higgs-induced attraction was discussed 
long ago for light quarks in in which case 
the critical number Nc for these objects to exist 
is astronomically large. In order to reduce this 
number, one should look for particles with a 
stronger Higgs coupling, i.e. a larger mass: thus 
this paper is focused on the heaviest members of 
the SM, the top-quark and the W-boson. The 
interest in the issue of "top bags" originated 
from the question whether a sufficiently heavy 
SM-type fermion should actually exist as a bag 
state, depleting the Higgs VEV around itself 
[S] . Although classically this seemed to be pos- 
sible, it was shown in refs (S1[TT1[T3] that quan- 
tum (one loop) effects destabilize such bags, ex- 
cept at so large coupling at which the Yukawa 
theory itself becomes apparently sick, with an 
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instability of its ground state. The issue rest 
dormant for some time till Nielsen and Froggatt 
[7] suggested to look at the first magic number, 
12 tops+antitops corresponding to the maximal 
occupancy of the lowest Z = 0, j = 1/2 orbital, 
with 3 colors and 2 from t + t. Using simple for- 
mulae from atomic physics these authors sug- 
gested that such system forms a deeply-bound 
state. In ref.[S] some of us have checked this 
claim and found that, unfortunately, this is not 
the case. While for a massless Higgs there are 
indeed weakly bound states of 12 tops, for a 
realistic Higgs mass they disappear. Further 
variational improvement of the binding condi- 
tions for the 12-quark system were discussed by 
Richard [51 who confirmed that 12 tops are un- 
bound for Higgs mass mn > Mc(12) ~ 50GeV 
or so. The issue has been also discussed fur- 
ther in the second paper of the same authors 
[TU] , in the direction of a hypothetical 4th gen- 
eration of SM superheavy fermions. In particu- 
lar, a possibility that baryons {N = 3) made of 
such fermions could be lighter than the sum of 
their masses, and possibly produced before the 
fermions themselves. 

Let us also mention a number of papers dis- 
cussing more general theoretical problems in 
the Standard Model. We wish to mention two 
of them, specifically tied to the top quark and 
its relatively large Yukawa coupling. These are 
the issues of (i) quantum corrections and (ii) 
the so-called "top condensation" . 

Regarding (i), the magnitude of quantum 
corrections to single-toj> bags has been ad- 
dressed by several authors. It is known that 
quantum corrections can deflate single-top bags 
into inexistence, or even turn the vacuum un- 
stable! (see, e.g. [TTtiTS] ). Similarly, the role 
of quantum corrections may be crucial to the 
existence of multi-qaanta bags. This led two 
of us [13] to study the magnitude of one-loop 
quantum corrections to such bags in a scalar 
approximation (which ignores the spin of the 
top quark and the Pauli principle) as a func- 
tion of N and g. Specifically, by writing the 
quantum-corrected energy of such states as a 
loop expansion 



E{N) = E,, 



l + ^Cfe(W)( 



Att' 



(1) 



we calculated the one-loop coefficient Ci(A^), 



following both the specific model and the prac- 
tical tools developed in 6 . The central is- 
sue was the A^-dependence for different cou- 
plings: explicit calculation of Ci{N) for g — 
1..8 showed that the for 5 ^ 1, the correction 
for large bags is dominated by the scattering 
phases of vacuum fields on the bag, and that 
Ci ~ 0(1) for large N without any strong N- 
dependence. For top quarks (the Yukawa cou- 
pling is g^/47r « .08) the one-loop quantum 
corrections are about IS.E — Q.OQEci.- We thus 
concluded that for Z, top-bags they are un- 
der good control, although for hypothetical 4-th 
generation fermions with a coupling a few times 
larger (i.e. mass > 1 TeV) quantum corrections 
get large. 

The second issue is the idea of top condensa- 
tion (for review and references see e.g. [H]). A 
(hypothetical) strong attraction in the tt chan- 
nel can lead to chiral symmetry breaking and a 
nonzero < tt > condensate. The lowest scalar 
meson would be identified with the Higgs bo- 
son. Whether the Higgs is composite or not 
is unimportant at low density, at the onset of 
multi-top bound states . However, as we in- 
crease the number of quanta and get deeper 
bound states the top density grows and one 
should confront the issue of the ultimate fate 
of the top bags. Two logically possible scenar- 
ios are (a) saturation at a certain finite pos- 
itive energy per particle, analogous to, e.g., 
nuclear matter, or (b) creation of a "chirally 
restored vacuum" with zero chiral condensate 
< tt >= 0. Fortunately, as we will show below, 
the W,top bags prefer configurations in which 
the density at large N is saturating, thus there 
seems to be no need to discuss a possible chiral 
transition and other issues related with it. 



C. The main approximation and the plan 
of the paper 

Trying to make the paper as simple as possi- 
ble, we ignore many secondary issues which can 
be easily included. In particular we ignore all 
forces other than Higgs-induced ones, such as 
e.g. strong interaction between top quarks or 
electromagnetic energy related to total charge 
of the bag. We set to zero the Weinberg angle, 
making Z degenerate with W. (In fact we will 
call below all gauge quanta W, for brevity.) As 
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we deal with large enough particle numbers N, 
one may assume that total electric, electroweak 
and color charges of the bag can be made suffi- 
ciently close to zero to ignore their mean fields. 
The lifetime of W, Z, t in the bag will not be dis- 
cussed in this paper, but it would be delegated 
to a companion paper [5]. 

The remainder of this paper is structured 
as follows. In section |lT] we give a qualitative 
discussion of multiple types and properties of 
bags. In section III we discuss W-boson bags, 
for which additional technical details are in- 
cluded in Appendi x \K\ Then, in section [TV] as 
well as in Appendix |B[ we study top quark bags 
and in section |V] bags containing both W's and 
tops. In section |VI| we discuss finite tempera- 
ture effects and we finally conclude with some 
comments about cosmological applications. 



Consider now the addition of fields with a 
global internal symmetry, leading to a con- 
served particle number N associated to each 
field. If these N particles (fermions or bosons), 
couple strongly to the Higgs field, or iV 1, 
the Higgs could be strongly distorted in the 
vicinity of these particles. To describe this 
situation, we adopt a mean-field approxima- 
tion, in which all the particles are described 
by the same wave functions in the background 
of a non-trivial Higgs field. Corrections to this 
mean-field description, such as, many-body, re- 
coil and retardation of the Higgs field are ex- 
pected to be suppressed by factors v/m, ran /m 
and 1/N. In the semiclassical approximation, 
the total energy of the system will be given by 
(see derivations in the Appendix A for PFand 
Appendix B for fermions) 



II. QUALITATIVE DISCUSSION OF 
BAGS 



A. Creating a bag in the Higgs vacuum 

According to the SM, the dynamics of the 
Higgs field $ is described by the Lagrangian 
density 



2 _ ^ 

4 ^ 



(2) 



where is the corresponding (gauge) covari- 
ant derivative. The usual vacuum state there- 
fore corresponds to a homogeneous Higgs field 
(<f>) ~ V permeating all of space. In the con- 
ventional unitary gauge, the Higgs field $ is 
represented by the real, dimensionless, field (j) 







(3) 



where v is the vacuum expectation value (VEV) 
which the Higgs achieves when (p — \. Assum- 
ing spherical symmetry, the energy reads 



E 



Higgs 



dr r 



A' 2 



(4)' 

where nijj = 2Xv^ is the Higgs physical mass. 
We will take v = 246 GeF and wih = 100 GeV 
and use units of 100 GeV throughout the paper. 



Eri — E 



Higgs 



E 



(5) 



where {Sa} is the spectrum of the correspond- 
ing field in the Higgs background, Ua is the oc- 
cupation number of each state and N — J2a 
is the total, conserved, particle number. 

In the Higgs vacuum, i.e. (j){r) = 1, the 
state of lowest energy with N particles corre- 
sponds to these particles sitting at the bottom 
of the continuum spectrum, with total energy 
NM. However, in the background of a non- 
trivial Higgs field there are two competing ef- 
fects. On the one hand, the gradient and poten- 
tial terms in ([5| increase the energy but, on the 
other hand, there might be some bound states 
levels with energy < Ea < M which can al- 
locate the quanta, lowering the energy of the 
system of particles at the expense of creating 
such distortion. Therefore, we consider ([s]) as 
a functional of (f)(r) and search for non-trivial 
bag solutions. 

Let us start by a crude estimate of the the or- 
der of magnitude of N for which such bags may 
exist. If we were to deplete a certain large vol- 
ume of the Higgs VEV (surface/kinetic terms 
neglected for now), it would require an energy 



Vol- 



m 



For a bag of radius, say, R ■ 100 GeV— 4, this 
energy is about 20 TeV. Thus, if the lowest 
W-boson energy level has a binding energy of 
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the order of 30 GeV per W or Z, an order of 
0(1000) of them would be needed to compen- 
sate for the bag energy and obtain some bind- 
ing. The top quarks are heavier and may get 
much larger binding, so one might naively think 
that less of them would suffice. But the situa- 
tion with fermions is more delicate due to the 
the Pauli exclusion principle: we will discuss 
this issue in detail below. 



B. Bosons in a bag 

The bags in which Higgs VEV is modified 
nearly to zero in some volume -to be called the 
"no-Higgs" bags below - makes particles which 
are heavy in vacuum nearly massless. Such 
bags are analogous in spirit to the various bags 
used in the 1970's to model the QCD confine- 
ment (e.g. the famous MIT bag [E]), except of 
course there is no need to take the mass out- 
side the bag to infinity. If the particles consid- 
ered are gauge bosons, W, Z, they all occupy 
the lowest bound state level. The energy of a 
bag with radius R is 



E : 



R 



+ CvR 



(6) 



where the first term is the kinetic energy of (ap- 
proximately) massless quanta confined in the 
bag. The precise constants depend on the de- 
tails of the model, one example being Cr « 
2.04 for a fermion confined in the MIT bag. 
It is clear that in such a case there is always 
a minimum and minimization with respect to 
R, at constant N, gives the size R{N) of the 
bag as a function of N. One finds that the 
total energy per particle decreases with N as 
E/N ~ iV~^/^, for particles with mass Mi, the 
binding {E/N < M) occurs above the critical 
number 



4 



^1/3 ^ 



(7) 



Note also that even if the bags are not ab- 
solutely bound, there might exist metastable 
bags, local minima with E/N > M, for TV* < 
N < Nc, while for N < N^ bosonic bags do not 
exist at all. 



There is another, maybe less obvious, possi- 
bility of the bag structure, which we shall re- 
fer to as the "inverted bag" . In this configura- 
tion, we let the Higgs reverse its sign inside the 
bag, crossing zero at some finite radius rg, and 
asymptotically taking the usual Higgs VEV v. 
If the Higgs field comes close to —v inside the 
sphere, such a bag is simply filled with another 
Higgs vacuum, so there is no volume energy but 
only a surface contribution. It is reasonable to 
expect (and we shall see later that it is indeed 
the case) that particles will be localized around 
the Higgs' node since they are massless there. 
If so, a large number of them would form a 
2-dimensional gas moving approximately on a 
sphere surrounding an opposite-sign vacuum. 



C. Fermions in a bag 

Let us start with a "no-Higgs" bag picture. 
As in any macroscopic situation, both the en- 
ergy of the gas and of the Higgs field scale as 
the volume and the (mechanical) equilibrium in 
such case is achieved when the internal pressure 
due to fermions is compensated by the exter- 
nal pressure of the Higgs vacuum. It is easy 
to see, on dimensional grounds, that for mass- 
less fermions the pressure scales as pf ~ /i^ 
and hence N = dE/d^ ~ ^^R^. Therefore, 
H ~ N^^^/R, and hence the total energy reads 



E = 



AvN^ 



CvR ■ 



(8) 



Minimizing the energy over i?, with fixed N , 
provides the pressure balance condition men- 
tioned above, leads to 



E _( c\/'Ay ' 
N \ 3 



3/4 



(9) 



The binding E/N < Mi, contrary to bosonic 
bags, gives a condition which is independent of 
N and whether binding will occur or not will 
depend exclusively on the values of the coupling 
constants and Higgs VEV, for any number of 
particles. 

Another option, an "inverted bag" , would 
put fermions at the Higgs VEV node at r « rp, 
within a shell of the width 0{l/mif). There- 
fore, for large N we expect the system to be 
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described by a 2-dimensional fermionic gas of 
volume cx Tq. For a 2D gas the pressure scales 
as p/ ^ fi^ and hence N ~ /i^^'o ^-iid therefore 
the total energy is given by 



with the associated conserved charge 



N 



d^xf. 



(13) 



E 



AsN^ 



Csrl 



(10) 



with two other constants As,Cs- The energy 
per particle 



E 
N 



= 3 



(11) 



is again saturated (N-independent). So 
whether binding occurs or not depends on the 
value of couplings and Higgs VEV. The pos- 
sible existence of the configuration may seem 
less intuitive, but a closer look at our sys- 
tem reveals an interesting situation, which sup- 
ports this picture for fermions more than for 
bosons. For very large bags of this kind, the 
system becomes effectively 1-dimensional (see 
Appendix [B]) and the Higgs equation ( [32| ad- 
mits a 1-dimensional kink solution, located at 
r — vq » l/rnu- It is well known that the 
Dirac spectrum in this background contains a 
series of discrete fermionic levels, including a 
fermionic zero mode, so in this limit the local- 
ization of (at least some) fermions on the (large 
sphere) surface costs no energy at all! 



III. GAUGE BOSON BAGS 

A. General equations for E, L and M 
modes 



The classical equations of motion read 

{a + M^cl)^)^ -d^d^W^^O . (14) 



By taking a derivative of Eq. (12), one can 



rewrite the equation of motion in the more 
transparent form 

(□ + M^^2)Ty^ + a^(^^^^^) =0. (15) 

It will prove convenient to study these equa- 
tions of motion in the electric (e), longitudinal 
(1) and magnetic (m) basis (see, e.g., [13]). As- 
suming stationary fields and spherical symme- 
try, we write 

^ie,un) ^ Yl:J^"'^f^,j^^){r)/r, (16) 

where Yj„i are spherical harmonic vectors and 
/e(r), fi{r) and fmif) are the radial wave func- 
tions for each mode. As seen in Appendix [A] 
in a static, spherically symmetric, background 
the last term in ( 15 ) vanishes for the mag- 



netic mode, leading to the simple Klein-Gordon 
equation for the radial wavefunction 



dr^ 



jU + 1) 



.f,n{r) = (17) 



for which j > 1- However, the Laplacian mixes 
the electro-longitudinal modes, leading to the 
set of coupled equations 



Consider the propagation of VF-bosons in an 
external Higgs field as described by the elec- 
troweak sector of the SM Lagrangian density 
(see Appendix [A]) 



d^ 
dr^ 



i{] + 1) 



/e(r)+(18) 



2y7(7TT)(l-r^)/Kr) 



= 0. 



£ = £ 



1 



Higgs 



Due to the remaining global symmetry Wi, — )■ 
e^"W^, there's an associated conserved current 
given by 



J' 



(12) 



2r 

dr 



jU + 1) 



M"^ ^ 



JjAr}dhi^\ 2^j{j + l)Ur) 
r dr 



/KK)i9) 



0. 



While these equations have similarities with 
the Klein-Gordon and Dirac equations (if the 
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latter one is written at the second-order dif- 
ferential equation), but there is also an impor- 
tant distinction. The last term in (15), which 
is not present for scalars or spinors, becomes 
singular in the region where the Higgs field </> 
crosses zero, effectively repelling the longitudi- 
nal modes from the node (or the whole bag, in 
the case of the "no-Higgs" bag) . It is qualita- 
tively explained by the fact that massless gauge 
fields have no longitudinal mode at all. As a 
consequence, below we will only discuss mag- 
netic modes, which are unaffected by mixing 
with the longitudinal one. 



f"- 

J m 



j{j + 1) 



and hence we find 



el=io'^ = Mw(j)'[ro) + 



3{i + 1) 



(21) 
(22) 



for the lowest level. Therefore, for large bags 
with a node, the wave functions of magnetic W- 
modes are Gaussian and localized around the 
Higgs's node, with energy eigenvalues given by 

m. 



B. Spectrum 



1. Analytical Results 



As mentioned before, the magnetic modes 
are not mixed with the longitudinal W boson 
modes which get repelled from the bags, so we 
study the spectrum of (17) delegated studies of 



the j — longitudinal mode in the Appendix. 
We shall separately study the cases of an in- 
verted and no-Higgs bag. 

In the case of a bag with no node, assum- 
ing that the Higgs takes some constant value 
(j)o with < </)o < 1 for r < i?, and <j) = 1 for 



r > R, the solutions to (17) are given by spher 



ical Bcssel functions with frequency squared 
cj^ — Af^0o- Assuming that the wavefunction 
is regular at the origin and vanishes at r = i?, 
the spectrum for the lowest radial excitation of 
the magnetic mode is given by 



1/2 



R 



(20) 



where /3 



is the I'st zero of the Bessel 



function and j > 1. Therefore, bound magnetic 
W-bosons will fill the whole of the bag, occu- 
pying the lowest level with an energy given by 
(1201). 



In the case of an inverted bag, we assume that 
4>(r) has a node at large tq, i.e. 0(^0) = 0, with 
0'('^o) > 0. We expand around x = r — — Q 
and Eq. (17) is a simple harmonic oscillator 



2. Numerical Results 

In order to find bag solutions for finite N , we 
adopted a variational approach and took as a 
trial function for the Higgs, the Gaussian profile 



") = 1 — a exp (— r^/w^), 



(23) 



where the two parameters, a and w describe 
the depth and the width of the bag potential, 
respectively. Solving the W-boson magnetic 



equation ( 17 1 in this Higgs background is rather 



straightforward using a shooting method, and 
we may compare the results to the analytical 
expressions for large-size bags derived previ- 
ously. For instance, for the case of an inverted 
bag (i.e. a > 1), the Higgs has a node at 
Tq = w^/\na and from (22) we find 



£0 




(24) 



for the lowest magnetic mode. In Fig. [T] we 
compare the numerical with this analytical re- 
sult for a bag with a = 1.3 and find good agree- 
ment. Similarly, we also found numerically the 
spectrum for a "no-Higgs" bag with a = 1. 



C. Existence of pure W-bags 

As shown in Appendix [Aj the total energy 
of the system with N W-bosons in the lowest 
energy level £0 is 



E 



Higgs 



Neo. 



(25) 
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FIG. 1: Lowest magnetic mode {j = 1) levels for a 
Gaussian bag with a — 1.3. The numerical results 
(points) are in good accordance with the analytical 
result 1241. 



Using a variational approach based on the 
Higgs ansatz (23), wc calculate the Higgs en- 
ergy (|4| in terms of a and w 



E 



Higgs 



(27raw)2 



H 



1 



3u; 
4 



303 



a 



8V2 



(26) 



and search for the local minima in this 2- 
dimensional parameter space. Since the behav- 
ior of bags with a > 1 and a < 1 is different, 
we present these cases separately. 



The no-Higgs Bag 



FIG. 2: Total energy/100 Gel/ of bag with a = 1 
(red curved line) and that of simple plane waves 
(black straight line) as a function of A'^. For A*' > 
Nc ~ 4700 the bags are stable, while for Nc > N > 
Nt ~ 1400, bags exist only as metastable states. 
For N < Nt: bags do not exist at all. 



plane waves (black straight line), as a function 
of TV. For iV > iVc ~ 4700 the bags are en- 
ergetically favored and thus stable, while for 
Nc > N > ^ 1400, bags exist only as 
metastable states. For N < N^, the bag would 
be to small to hold any levels and cannot exist 
at all. As one can see from this figure, the gen- 
eral feature of the bags in question is that the 
overall energy per particle is not very different 
from the particle mass: and yet the way this is 
reached is completely different from individual 
particles without a bag, since the quanta are 
near-massless inside the bag. 



For our cosmological applications it would be 
of particular importance to have near-empty 
bags with aw 1, because only there the 
electroweak sphalerons can be represented by 
pure-gauge (COS) sphalerons which have larger 
probability. So, in order to demonstrate this 
situation, we took a — 1, solved the W-boson 
spectrum numerically and proceeded to calcu- 
late the total energy of a bag with N quanta. 
In Fig. [2] we show the total energy of bag with 
a — 1 (red curved line) and that of simple 



2. Inverted 



From Eqs. 
inverted-type 



(24 



26) we also searched for 
bag solutions by look- 



[a > 1) 

ing for local minima of the total energy in both 
the a and w directions at fixed N . The bind- 
ing of such bags occurs for a rather larger N, 
but do occur. In Fig. [3] we show an example 
for N = 20, 000 for which a local minimum was 
found at a = 1.4 and w = 7.5. 
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FIG. 3: Total energy/100 Ge\/ of an inverted-type 
W bag with A'' = 20, 000 in the lowest magnetic 
mode in terms as a function of a and w ■ 100 GeV. 
In this example, a local minimum in both directions 
is found for a = 1.4 and w — 7.5. 



IV. FERMION BAGS 

A. The Dirac equation 

Omitting for now gauge bosons, we consider 
a system of N heavy fermions interacting with 
a background Higgs field as described by the 
SM Lagrangian density 



1 / F(r)r!,,„ 

r [ (-l)i/2(i+'-'')G(r)%p. 



(29) 



where ^jim are spherical 2-component spinors 
and we take normalization J dr [F"^ + G^) — 1. 
The so-called Dirac parameter k, is defined as 



I 



for j = 
for j 



1 + 1/2 
-1-1/2 



(30) 



and runs over all nonzero integers, being posi- 
tive for anti-parallel spin and negative for paral- 
lel spin. Dirac's equation reads (See Appendix 
iBl): 



[e — m<j)) F 
(e + m (/)) G 



-G' + {K/r)G (31) 
F' + [n/r) F 



The form of these equations presumes that the 
eigenvalue e is positive. A negative eigenvalue 
would correspond to a state in the lower fermion 
continuum. If so, a charge conjugation trans- 
formation turns it into a positive eigenvalue for 
an antifermion. The Higgs equation of motion 
reads 



' + ^0(1- 



{N - l)m F^ - G^ 



An v'^ 



(32) 

Note that in the r — ^ cx) limit, the source term 
in the right hand side, as well as the additional 
4>' /r term from the Laplacian can be neglected 
and the equation becomes the usual equation 
for a ID kink. 



C = Cmggs + ^ (^7^9^ - 5*) ^, (27) 

where is complex scalar in the fundamen- 
tal representation of SU (2) ^ and ip is a, Dirac 
spinor in 3-1-1 dimensions. Due to the global 
symmetry ip e*"?/", there's an associated con- 
served charge given by 



N = J Sxtphdoil), (28) 

which is identified with the fermion number. 
For completeness, let us remind the reader the 
standard notations for Dirac spinors in spheri- 
cal coordinates [T7] : 



B. Spectrum 

1. Analytical Results 

In order to analytically justify the physical 
picture of inverted bags , we calculate the en- 
ergy of surface-fermions for an inverted bag of 
large (but finite) radius. We assume that (j){r) 
has a node at large rp, at which the derivative 
is positive, 4''{ro) > 0, we search for solutions of 
the Dirac equation in the near-Gaussian form 

F{r) = A(r)e-^M ^ (33) 
G{r) ^ S(r)e-^M , 
S{r) = m r(j){r')dr' ~ ^ m (/)'(ro) (r - tq)^ . 
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The large mass m here favors the locahza- 
tion of F{r) and G{r) in the vicinity of rg, 
which justifies the approximation. The func- 
tions A{r), B{r) can be iteratively found from 



eA- 
eB 



B' - {K/r)B 
-A'-{K/r)A 



mcj) {A - 
-m(j) (A - 



B) 
B) 



(34) 



Since A, B arc smooth functions of r we expand 
to first order about r — rg and derive a system 
of four hnear homogeneous algebraic equations 
for (^0, Bq, A'q, B'g), given by the matrix 



Dirac spinor), which has a particular chirality 
(i.e. F = ±G) and thus no scalar source. How- 
ever, this is only true for the leading terms: As 
To becomes smaller, the ^ l/j'o terms of our 
expansion violate this equality, rendering the 
scalar — nonzero. 



2. Numerical results 



K 







ro 




£ 


-1 


'A 


m<f)'{ro) + 4 
'o 


e 


K 


m(f)'{ro) 


K 

ro 


^ J 



Setting det — 0, we solve for e and get 



= m0'(ro) + k'^/tI ± . m'^(t)''^{ro) + 



1 

74- 



V '0 

(35) 

We see from this expression that the levels e+ 
are finite for r — )■ oo, while e_ goes to zero. The 
latter correspond to the would-be zero modes 
in the 1 dimensional case. Therefore, for large 
enough bags these will be below any other level 
(see Appendix |b] for explicit results) and in 
what follows we consider these levels only. Note 
that to order O(^), 



''0 



(36) 



hence in this approximation the spectrum is in- 
dependent of the mass of the fermion and the 
value for the first level {k=-1) is less than 2.04, 
which appears for the volume bag. Thus, not 
only the Higgs energy, but the fermionic energy 
benefits as well. Note also that to this order, 
the functions A^ B are simply constants whose 
values are fixed by the normalization condition 
to 



A 



B = [m(/)'(ro)/47r]i/'* 



(37) 



Therefore, in this limit F^ ^ and thus these 
states do not contribute to the scalar source 
F^ — G^ in the Higgs equation. This is be- 
cause this approximation corresponds to the 1- 
dimensional zero-mode solution {F and G be- 
ing the upper and lower components of a 1+1 



In this section we provide some of the numer- 
ical results, and comparison to the analytical 
results derived before, for the the Dirac spec- 
trum in a Higgs background. We took an ansatz 
for the Higgs which interpolates between the 
type of bags discussed above, namely surface 
and volume bags, parametrized as 



(r) — {1 — rj) + rj tanh 



Ar 



(38) 



where < 77 < 1 and < i? < 00, with 
7] = corresponding to the trivial Higgs vac- 
uum, ?7 = 1/2 to a no-Higgs bag and 77 = 1 to 
an inverted bag. By using a shooting method, 
the spectrum for Dirac's equation is found nu- 
merically and examples of the levels for an in- 
verted bag with r] — 1 are shown in Fig. [4j 
This bag has a node at rg = R. One can see 
that there is very good agreement between the 
analytical expression ( 35 1 and the numerical re- 



sults (points) . Although it has been derived for 
large R, the agreement is not there only for 
small R < 2 which we will not discuss in this 
work. 

The Table shows magic numbers and the or- 
der in which levels are populated. Some levels 
are also shown in Fig. [s] we show some levels 
(up tol = 2). As the bag's size decreases (or for 
smaller 77), some levels with k > 0, and radial 
excitations of k < states, cross and one must 
carefully fill the levels accordingly. Note that in 
this case, the lowest levels with k < asymp- 
totically become zero modes (with the correct 
chirality), while radial excitations thereof and 
levels with k > merge, asymptotically becom- 
ing the next fermionic discrete level. 
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FIG. 4: Energy levels (in units of 100 GeV) of the 
Dirac operator in the background of an inverted bag 
(j) = 1) as a function of the bags's size R - 100 GeV . 
These levels correspond to increasing I (from the 
bottom up) for k < only. The analytical e_ re- 
sults for large R (solid lines) are in good accordance 
with the numerical results (points) for even rather 
small R. 



FIG. 5: (color online) Dependence of some bound 
state levels energy e in units of 100 GeV on the 
size of the bag, expressed as the parameter R in 
units 1/100 GeV~^ for a bag with r] = 1. The color 
coding of the levels and their quantum numbers are 
listed in Table I. 



C. Non-existence of pure top bags in 
vacuum 
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Deg.{tt) color 


1 





-1 1/2 


12 


blue 


2 





-2 1 3/2 


24 


red 


3 





-3 2 5/2 


36 


green 


4 





1 1 1/2 


12 


black 


5 


1 


-1 1/2 


12 


blue 


7 





2 2 3/2 


24 


yellow 
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1 


-2 2 3/2 


24 


red 


9 





3 3 5/2 


36 


violet 


10 


1 


-3 2 5/2 


36 


green 



TABLE I: The properties of some levels, including 
the number of radial nodes rir, Dirac parameter 
ft, orbital momentum I, total angular momentum 
j, multiplicity of states Degiior tt bags) and color 
code used in our figures. Obviously, for pure top 
bags the multiplicity is halved. 



1. No-Higgs Bag 

Consider as a trial configuration a spherical 
cavity of radius R with constant zero Higgs 
VEV inside, a small transition region of size 
l/niH, and a constant Higgs VEV outside, just 
as in the Friedberg et al. bags, except that we 
shall occupy not only the lowest level. Note 



that this is a solution to (32 1, except in the 
(small) transition region. This will correspond 
to a free gas of massless fermions in a spherical 
volume of radius R and statistical mechanics 
tells us that 



= {6n/D)'/^ (3/4) 



^4/3 



R 



(39) 



with D the degeneracy factor due to spin, color 
and particle/antiparticle. From here we read 
the coefficient Ay in ([8|. Directly from ([S]) we 
also find Cv = ™w and therefore the en- 
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ergy per particle reads 



E 
N 



D 



1/4 



(40) 



The bag will be bound provided E/N < rrit- 
Clearly, having a larger number of fermionic 
species helps: say t + t bags, with D = 12, are 
more economic than pure top bags with D = Q. 
From this expression (assuming v = 246 GeV) 
we find that the bag will be bound provided 
that ruH < 77 GeV which, unfortunately, is not 
the case in vacuum. 



2. Inverted Bag 

Given the scmiclassical solution derived in 
|IVB H we can estimate the binding energy for 
very large inverted bags. Finite-size bags re- 
quire a detailed study of all the fermionic lev- 
els, including radial excitations and k > 
states. However, as seen in Appendix |B] for fi- 
nite but large enough bags the lowest levels are 
those given by (361 and therefore we shall fill 



these levels only for now. The maximum occu- 
pancy of each level is given by = D\k\, with 
D = 6, 12 for only t or tt bags, respectively. 
Therefore, for a large size bag the number of 
levels that it can hold, Kmax ~ fn^o, is large 
and from (|36| the fermionic energy reads 



E,i 



D (2iV)3/2 
3 ro 



(41) 



where we have used — '«'^/3 and N — 

Z?^|k| ~ Dk? /2. Thus, in the notations 
used for our generic estimates of the "surface 
bag" before, we have explicitly justified the 2-d 
Fermi-gas shell picture and, furthermore, fo und 
As = 23/2/3. For large r, equation ([32| 

is solved by the kink profile 



(r) = tanh[(r - ro)m/^/2]. 



(42) 



which leads to Cs — 8.38mHV^ . These values 



for C5 in ( 11 1 lead to an energy per-particle 
of 



E 
N 



3.68 



1/3 



(43) 



Again, we see that binding in principle is 
possible provided the values of the parameters 
are such that E/N < rrit. However, this 
requires the even more restricting upper bound 
rriH < 20 GeV, which is excluded experimen- 
tally and hence no binding of inverted top bags 



These results indicate that, for realistic Higgs 
mass, pure top bags are not energetically fa- 
vored over free tops. However, one may re- 
lax this condition and still wonder if at least 
metastahle top bags exist. To answer this we 
adopted a variational approach using the trial 
profile ( 38 1 and searched for local minima in 
(77, R). In the limit R ^ Ar and using Q, 



we write the total energy in terms of rj, R and 
AR. From here we find that the optimal value 



for Ar is Ar(°P*) = 



mH(2~ri) 



which leads to the 



expression 
27rw2 



R^ruHv'' (2 - v) 
+ ^R^mW (1 - vf 



(44) 



For r] = {0, 1}, corresponding to the usual 
Higgs vacuum and inverted bag, the volume 
term vanishes. If 77 > 1/2 the bag has a node 
at ro = i? -I- (2/mij) tan~-^[(77 — l)/r]] and from 



from ( 35 1 we can parametrize the fermionic en- 



ergy in terms of (77, R) as well, which allows us 
to search for local minima in the 2-dimensional 
parameter space. Although an interesting spec- 
troscopy of these objects are found for small 
Higgs masses (e.g. tuh ~ 50 GeV), for a real- 
istic Higgs mass, no local minima were found 
for any TV. Therefore, no pure top bag solu- 
tions, stable or metastable, were found in our 
calculations. 



V. TOP QUARKS IN W-BAGS 

Although we have seen that pure top quarks 
are not heavy enough to create top bags by 
themselves (in vacuum), W-bags containing 
N > Nc W-bosons certainly do exist and given 
that top quarks are present, one should address 
how is it that fermions behave in the presence 
of these bags. From the results of the previ- 
ous section, we know that fermions will bind to 
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these bags, forming a 2-d fcrmi gas in the sur- 
face of an inverted bag or, fiUing the volume of 
a no-Higgs bag. The next issue to study is the 
relative position of the W and top levels. 

Free (or weakly bound) top quarks are much 
heavier than W bosons and thus decay into an- 
other quark and W. In the presence of a large 
Higgs bag, however, the relative position of W 
and top levels can be such that quark levels 
are lower than gauge boson ones. If this is the 
case, the life-time of i-quarks in such bags will 
be much longer, given by next order decays into 
three fermions, like in the usual beta decays. 

To answer when this happens we numerically 
studied the behavior of these levels as a func- 
tion of the bag's size. In Fig. |6] we show 
the resulting energy levels for the M-mode (red 
points/line) and the top quark levels (small 
black circles) in the background of a Higgs bag 
with a = 1.2. Indeed, it is observed that when 
the bag becomes large the lowest-lying quark 
levels are below the lowest W ones, preventing 
their decay in the lowest order of weak interac- 
tions (to Wb). We discuss the possible cosmo- 
logical role of these observations elsewhere 0] . 
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FIG. 6: The energy of the levels Ew/ WO GeV ver- 
sus the size of the bag w- 100 GeV , for the Gaussian 
ansatz for the Higgs and a = 1.2 (with zero case). 
Black (small) circles are some fermionic levels with 
K < and the large (red) circles correspond to the 
W magnetic level for j = 1. 



VI. FINITE TEMPERATURE 

As it has been already emphasized in the In- 
troduction, these multi-quanta bags can hardly 
be produced experimentally, and their main ap- 
plication is cosmological, as "doorway states" 
facilitating electroweak sphaleron rates as well 
as possible CP violation, hopefully to result in 
a possible mechanism for baryogenesis. In gen- 
eral, the possible scenarios can be divided into 
(i) Standard Cosmology, with very small devia- 
tions from thermal equilibrium and (ii) hybrid 
models, in which the electroweak phenomena 
are coincident with the end of inflation. 

Although our paper [3] is devoted to the lat- 
ter scenario, with strong deviations from equi- 
librium, in this section we will restrict ourselves 
to the thermal equilibrium case, as a necessary 
intermediate step, extensively studied in the lit- 
erature. For an overview and references see, 

e.g., m- 

Before we come to specifics, let us outline 
the main issues to be discussed in this section. 
As the temperature in the early universe is "at 



the electroweak scale", with momenta of the 
order of 100 GeV or so, all parameters of the 
SM are significantly renormalized. This in par- 
ticular refers to the Higgs mass Mh{T) and its 
VEV w(r), which are strongly reduced as T ap- 
proaches the critical temperature Tc- 

Recall that, as emphasized many times 
above, if the Higgs mass would be smaller, 
it would significantly increase binding of the 
multi-quark bags. Therefore, one may naively 
think that a high temperature environment will 
foster the existence of W-top bags. However, 
the part of the top and W masses related to 
Higgs are also reduced, and so is their binding 
to the bag. On the other hand, only a part of 
the particles's masses come from Higgs mech- 
anism: at finite T they are complemented by 
the so called "thermal masses" due to particle 
rescatterings. The outcome of these two effects 
is the subject of this section. 
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A. SM parameters 

A detailed review and discussion of the gen- 
eral issue of finite-T electroweak theory and its 
"phase transition" from symmetric to the bro- 
ken phase is obviously beyond the aims of the 
present paper. Let us just say that even in 
the relatively weakly coupled electroweak the- 
ory one still finds elements of a strongly inter- 
acting gauge magnetic sector, and other non- 
perturbative issues, which are reminiscent of 
QCD phase transitions. Here we limit ourselves 
to the discussion of much simpler effects, such 
as the lowest-order rescatterings or "thermal 
masses" , which are dominant, as we are not too 
close to Tc- We will begin with the Higgs boson 
and then proceed to gauge quanta and finally 
to fermions. 

The temperature dependence of the effective 
Higgs potential Vef f{v,T) is quite an involved 
subject, which has been studied by renormal- 
ization group methods as well as numerical 
lattice simulations. Let us just say that at 
high temperatures the Higgs mass and VEV 
are significantly renormalized, eventually melt- 
ing away at a critical temperature T = Tc. For 
practical applications, not too close to Tc, the 
Higgs mass and VEV scale down in the same 
way, approximately as 



= 1 - 



y2 



(45) 



The masses of W-bosons and quarks receive ad- 
ditional contributions to which we turn now. 
The first effects are lowest order rescattering 
given by the single one-loop diagrams. For 
gauge bosons at T higher than all masses the 
result is [TS| 



(46) 

where the first term is the usual vacuum mass, 
but with a T-dependent VEV, and the last 
three terms come from W rescatterings on 
gauge, Higgs and (the single species) quark, by 
which we mean the top quark. In fact these are 
due to lowest-order tree scattering diagrams. 
Substituting v{T) from (45) one finds that 
the negative coefficient from a decreasing v{T) 
and positive coefficients from rescatterings 
tend to cancel each other! (In fact they do 



when Tc w UQGeV.) Thus we wiU adopt 
an approximation in which Mw, unlike toj^, 
does not have a significant dependence on T, 
due to the (diminishing) Higgs part nearly 
compensated by the (increasing) "thermal" 
part. We will address how this will affect the 
binding of W^s in the next section. 

For (top) quarks the description of their ef- 
fective mass is more complicated for several 
reasons. Similarly to finite-T QCD, the ther- 
mal masses due to rescattering on gauge bosons 
lead to chirally diagonal T,ll.rr contributions, 
while their interaction with Higgs generates the 
usual left-right LR mass term M . The result- 
ing Dirac operator in the chiral basis has the 
structure 



D = 



M+ S 



(47) 



RR 



with the "masses" with different chiral struc- 
ture. In contrast to QCD, weak interactions 
produce the LL term only, while strong leads 
to symmetric LL -\- RR contribution. Further- 
more, for quasiparticles with nonzero momen- 
tum p one has additional splitting in helicity 
cr = (J • p = ±1 into "quark" and "plasmino" 
modes as discovered by Klimov and Weldon in 
1970's, leading in total to four distinct modes 
per quark fiavor. For details the reader may 
consult e.g. the paper by Farrar and Shaposh- 
nikov |19j . from which we took specific expres- 
sions to be used. At small momenta splitting in 
helicity goes to zero and there remains two val- 
ues of masses, split in chirality, ^ll,rr(j-^,'P = 
0) = -^Il.rr/^^ where 



IrpZ 



27ra, 



(48) 



^hR —T ^ + —T , (49) 



while the LR mass term is still given by the 
Higgs by 



(50) 
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for generic up U quark flavors. If is a func- 
tion of coordinates, representing the bag, these 
expressions provide the form of the Dirac eqn 
for the bound states. 

Now, in the absence of gluons, and specifi- 
cahy for the top quark one finds that the Higgs 
mass part is much larger than the scattering 
part 



2 

LL 



O.IT^ 



(51) 

which in turn is larger than the RR compo- 
nent. Therefore, the coefficient of coming 
from the Higgs term is factor 20 larger than 
that coming from rescattcring. If the latter is 
neglected, one finds a simple approximation for 
the tops, in which all T-dependence of the top 
quark comes from the Higgs sector. If so, se- 
lecting e.g. v{T) to define a scale in a binding 
problem, one finds Mh{T), Mt{T) to (approx- 
imately) scale in the same way: which means 
that the Dirac equation for the bound levels re- 
mains to a good accuracy unchanged, except 
for the absolute scale. Thus all the results we 
reported above about tops in the bags remains 
valid, among them the unfortunate conclusions 
of absence of pure top bags, at any N, and thus 
they must coexist with certain number of the 
gauge quanta to form a bag. 

Note that in the previous discussion we have 
considered the gauge and Yukawa coupling con- 
stant to be constant: for completeness, let us 
comment now on what is known about their 
running. The running of the effective coupling 
in the 3-d magnetic sector at small momenta k. 
Its renormalization group flow is determined by 

EQl 



dln{k) 



24^^^'^'^) -k 



(52) 



which gets strong at momenta smaller than 
some critical value k < kc, where the 3d mag- 
netic theory goes into a non-perturbative con- 
fining regime. However this is only become im- 
portant very close to Tc and thus will not be 
included. 

The next issue is the T-dependence of the 
Yukawa couplings, especially the largest one 
defining the top quark mass gt- As noticed 
by Marciano [H], its beta function has two 



terms with an opposite sign: positive from self- 
coupling and negative, containing a^, from a 
virtual gluon. They tend to cancel each other, 
producing a quasi-fixed point. Furthermore, as 
a coincidence, the fix point value is not far nu- 
merically from the physical value of the cou- 
pling for the top quark (but not other quarks). 
As a result, one can neglect the running of the 
top Yukawa coupling. 

Our discussion in this section so far has been 
limited to expressions derived in thermal equi- 
librium. However the specific application we 
will discuss in |3| corresponds to the so called 
hybrid scenario in which there are strong de- 
viations from equilibrium: thus we would like 
to comment how those can be included. Some 
of the non-equilibrium effects - for example the 
fact that at the time under consideration gluons 
are not yet effectively excited - can be taken 
care of by simply omitting the corresponding 
contributions. Some others can be included 
by treating differently the "bulk" and the "in- 
side" of the WZ — topbags, by using the same 
thermal expressions but with different temper- 
atures and densities. In the approximation in 
which the scattering amplitudes can be consid- 
ered to have a weak enough dependence on the 
momenta, so that they can be factored out of 
the thermodynamical integrals, one can simply 
use the hard-thermal loop expressions, with 
substituted by by bosonic or fermionic integrals 



rj-iZ 



24 



2E (27r) 



3 '"P 



(53) 



I f ermionic 



48 



d^p 



2E{2tt 



(54) 



in which the occupation number Up is 
Bose/Fermi distribution including particle 
masses and non-zero chemical potentials. In 
particular, one can treat this way the top-gauge 
scattering inside the bags, in which the density 
of tops and W, Z quanta are created occasion- 
ally, by large fluctuations in the Higgs field in 
process of its primordial formation. 
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B. W bags at finite T 

As we already discussed in the preceding sec- 
tion, at finite T there are two terms in the ef- 
fective gauge boson mass, 

M^ = ^v^{T)cl){r) + M2iTf 

where the first term is Higgs-induced and the 
second term is generated by rescattering on the 
medium. The first term is due to the Higgs cou- 
phng and therefore space-dependent via the bag 
profile function (/'(r). The second term is how- 
ever assumed to be just a position-independent 
constant, which therefore can be absorbed into 
the effective frequency 

uj^ ^ = Lu^ - AhiT)^ (55) 

in all the equations of the preceding section. 
Furthermore, by using v(T) as the unit of scale, 
we are left with universal equations which have 
only T-independent profiles, from which uni- 
versal binding for Lb/v{T) combinations follow. 
However, while minimizing the bag, one has to 
return back to the usual energies. 

We noticed before that in the high-T 
quadratic approximations, the 0(T^) terms in 
the first and the second term tend to approxi- 
mately cancel each other. However the rescal- 
ing we are discussing now works for any T- 
dependence, and we do not use this specific ap- 
proximation. The expression for the W mass 
we use (in units of — (100 GeV^)^ as else- 
where) is ml^{T) = .24T2-I-.64, so the absolute 
magnitude of the W mass is growing with T. 
Regarding the T-dependent decreasing Higgs 
mass, this effect is even more pronounced. 

Using this approximation one can convert the 
universal spectrum of u) for the lowest j=l mag- 
netic modes in a bag into absolute units. FigjT] 
shows a sample of the results. As one can see, 
the binding energy is reduced as the tempera- 
ture grows. In our cosmological application in 
|3] we will have T^uik ^ 50 GeV, corresponding 
to the second curve from the bottom, which is 
to be compared with the zero T case shown by 
the bottom curve. 

Substituting these results into the expression 
for the total energy, together with the Higgs 
bag energy, one can look for a minimum over all 
parameters of the ansatz used, as before. What 



we found is that the change of the W mass does 
not significantly change the pattern of W bag 
binding. As N increases, W-bags develop first 
a local minimum for some N > N^, and become 
bound for TV > iV^ ~ 0(1000). 

A more important problem for the stability 
of the bags is of course "self-ionization" , i.e., 
the thermal excitation of a particle from the 
lowest bound level to the continuum spectrum. 
For low temperatures, such process will be sup- 
pressed by a Boltzmann factor. As T increases, 
however, the binding energy becomes compara- 
ble to the temperature and self-ionization is no 
longer suppressed, becoming quite substantial. 
Indeed, for T ~ 50 GeV (the second curve from 
the bottom in Fig[7]) one finds that the bind- 
ing energy per W is ^ 30 GeV. Therefore, the 
prospect of having long-lived W — Z — t bags 
at higher T is slim. This is why we focus on 
"hybrid" or "cold" scenario of baryogenesis in 




4 6 8 10 12 14 16 18 20 



FIG. 7: The lowest j = 1 magnetic modes in a 
Gaussian bags versus its size w. The four curves, 
from bottom up, correspond to T/(100GeV) — 
0, 0.5, 0.75, 1, respectively. 



VII. SUMMARY AND DISCUSSION 

We have studied the existence of multi- 
quanta bags in the electroweak sector of the 
SM. We found that pure W-bags exist provided 
the number of gauge bosons is greater than 
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Nc ^ 0(1000), with typical sizes of the order of 
R ^ O.lGeV^^^ . They can be either of the no- 
Higgs type, with W-bosons occupying a volume 
where the Higgs VEV has been depleted from, 
or of the inverted-type, with W-bosons local- 
ized in a surface shell of width 0{l/mH) where 
the Higgs vanishes, enclosing an inverted-sign 
Higgs vacuum. Although pure top bags are ex- 
cluded from the vacuum, top quarks will cer- 
tainly bind to these W bags and we found their 
energy levels in such bags. Envisioning cosmo- 
logical applications in our companion paper ^ , 
we have also studied such objects in a finite 
temperature environment. Although the main 
properties of such bags are unchanged at finite 
temperature, as temperature grows the ioniza- 
tion of these objects becomes important, even- 
tually evaporating them into disappearance at 
temperatures close the electroweak scale. 
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Appendix A: VP'-bosons in a Higgs 
Background 



transforms as a doublet under SU{2). Taking 
the unitary gauge one can present it via one 
real component 







(A2) 



Assuming that the scalar field develops the vac- 
uum expectation value = 1 and the Higgs 
mechanism takes place, one finds that the gauge 
field can be presented as a new C/(l) field A^, 
and a triplet of massive fields , 

Af, = -smeAl + cos9Bf,, (A3) 
= cos OAl + sin 9 , (A4) 
W, = (4 - ^Al) /V2 . (A5) 



Here = represents the ly-boson with 
electric charge e = — |e|, and W* 
boson with charge e 



„ W+ the W- 
II 

e and 6 is the Wein- 



berg angle. Expanding the Lagrangian Eq.( Al I 



in the vicinity of Z = A = and retaining only 
bilinears in the fields W^.W'^^ terms, one de- 
rives an effective Lagrangian 



1 



+M^(l3^w;w'' 



3+1 dimensions 



which describes the propagation of M^-bosons 
in an external Higgs field. As mentioned in the 
text, due to the global symmetry W,^ — )■ e*"VF^ 
one finds the conserved current 



a. Lagrangian, Conserved Charge and Equations 
of motion 

Consider boson fields in the electroweak sec- 
tor of the Standard Model (see, e.g., Ref.[lll) 

£= - J {d^A, ~ d.Af, + g X A,)\{Al) 

Here A^ and B^ are the triplet of SU{2) 
and the U{1) gauge potentials respectively 
(abridged notation is used here). The covari- 
ant derivative D^^^ takes into account that the 
Higgs field $ has a hypercharge Y = 2, which 
describes its interaction with the U{1) field, and 



J' 



Assuming that the fields are stationary with 
frequency uj, i.e. 

W,{x,t) = e-'^*W,ix) (A6) 

one finds explicitly the charge density 

f = -2uj\W,\^ -i{W'd,W^ -W*d'Wi)p) 

From the Lagrangian one derives the classical 
equation of motion for gauge bosons 

{a + M^cl)^)^ -d''d''W^ = . (A8) 
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Taking a covariant derivative in Eq.(14) one 
finds 

d^^W) = (f'd^W^' + = .(A9) 



Evaluating d^W^ fi-om Eq.(A9l and substitut- 
ing the result back into Eq. one rewrites the 
latter one in a more transparent form 



0. (AlO) 



c. Spherical coordinates 

We consider a static, spherically symmetric 
Higgs field <j){r). The angular dependence of a 
vector particle wave function will then be given 
by the spherical vectors (see, e.g., 



.(e) 



(A16) 



r{m} 
jm 



n X Y, 



■jm 



This form of the equation of motion will be 
more suitable for studying the behavior of dif- 
ferent W modes close to zero's of 6. 



b. Hamiltonian 

The Hamiltonian density is given as usual by 
the Legendre transformation 



(All) 



where pi, p* are conjugate momenta and are 
given by 

= -d^W*+d^W^* (A12) 
p* = -d^Wi + d,W°. (A13) 

and Pq = Pq = Q as expected. Using the equa- 
tions of motion ( A8 1 and performing some in- 



tegrations by parts, one finds 



uj 2uj\Wi\'^ + i{W''diW^ (AM) 
W"d,W°) 



Comparing to (A7| one finds Hw 



Fi- 



nally, integrating over space 



(A15) 



which is the desired relation. In order to eval- 
uate the total energy of the bag we must find 
the eigenvalues lj from the equations of motion 
for the occupied modes and also add the pure 
Higgs terms in the Hamiltonian (the kinetic en- 
ergy plus V{(j)) which has not been present in 
the formulae above. 



Here F, 



Yjfn{9,ip) is the spherical function 

and, Y^-^jYj-^^^Yj^'' are the electric, longitu- 
dinal and magnetic vectors. The symbol c?„ in 
Eq.(A16l indicates the angular part of the gra- 
dient, dF{e,(j)) = dnF{e,(j))/r, and n = r/r is 
a unit vector along the radius vector. Relevant 
formulas for the Laplace operator read 



"m 



]{3 + l)Y:l+'i^/W+l)Y^ 



(0 



A„l^;„V = 2v/7(JTT) - + 1) + 2)y^ 



(m) 



(m) 
jm 



(0 
jm ' 



Here A„ describes the angular part of the 
Laplacian, i.e. AF{9,(t)) = A„F/r2. The par- 
ity for electric and longitudinal polarizations 
equals P =^ (^l)''; for magnetic polarization 
the parity is P = (—1)^+^. The orbital mo- 
ment I takes the value I = j for the magnetic 
polarization, in agreement with the parity for 
this state (states pi, ^2, /s ■ • •)• 

The electric and longitudinal polarizations 
are constructed as linear combinations of the 
two states with Z = j ± 1: this is why they are 
mixed by the Laplacian operator. For j = 
there exists only one spherical vector, which is 
purely longitudinal and has 1 = 1. The mixed 
electric-longitudinal modes correspond to the 
following states with j > 1: Si,p2, ^i, ^3, /2 • • ■• 
The states with j = are given by the longitu- 
dinal mode alone (states po). 

For a static, spherical Higgs field 0(r), Eq. 



(A9| gives 



9.iy^ = 2Wn^. 

ar 



(A17) 



The electric Y- 



and magnetic Yj™'' vectors 



are orthogonal to n (Wn — 
longitudinal wave Wn = W. 



0), while for the 
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0.6 



0.4 



y=0,Z 



FIG. 8: The energy of the levels e/(100 GeV) versus 
the size of the bag w ■ (100 GeV\ for the Gaussian 
ansatz with a = 0.9 (no zero case). The j — 1 
M-mode is below the j = L mode. 



e. Electro-longitudinal modes 

For the longitudinal mode the effective po- 
tential ) remains and becomes sin- 
gular, ^ 1/0^, near cf) — 0. It is easy to un- 
derstand the reason for this singular behavior. 
The longitudinal wave does no exist for a zero 
mass particle. Therefore, the longitudinal wave 
must vanish in the area (f> = since the Higgs 
field (j> plays a role of the effective mass. 

It is convenient to use the substitution 



L(r) 



ipjr) 
r(j){r) 



(A18) 



to recast Eq. 
Schreodinger-like form 



( 19 1 into the following 



y(r)V'(r) 



+ 2<P{r)^j{j + l)Ur)/r') (A19) 



d. Magnetic modes 



For a magnetic polarization, Eq. (A17) tells 
us that dfj,W^ = 0. Therefore, the last term 



in Eq.(15) vanishes and fm satisfies the simple 
Klein-Gordon equation (17) with j > 1. Al- 



though also true for the electric polarization, 
due to the fact that the Laplacian mixes the 
electric and longitudinal polarizations one finds 
the set of coupled equations (18 19). As men- 



tioned in the text, this fact makes us expect 
that the magnetic modes are below the electro- 
longitudinal ones. In order to verify this, we nu- 
merically solved the W-boson spectrum in the 
background ( [23| for various values of a. The re- 
sults for a particular example with a = 0.9 are 
shown in Fig. |8] Comparing the j — 1 mag- 
netic mode (large circles) with the j = lon- 
gitudinal mode (large boxes) one can see that 
the former levels are lower, in spite of having 
a larger j. This is the effect of the repulsive 
potential, which moves the longitudinal states 
upward, even when no zero and singularities 
are present. Therefore, we fill only the lowest 
(j — 1) magnetic level when searching for W- 
bag solutions. 



with the last mixing term and the effective po- 
tential 



V{r) = 



^ h0(r) M^- 



1 



dr^ (j){r) 
{d(j){r)/dr)f 



dr 



(r) 



It now includes a new term with the first deriva- 
tive of the kink: yet the most singular term at 
the position of a zero is with the square of the 
first derivative of 4>{r). This l/(r — tq)^ sin- 
gularity makes the wave function vanish at tq, 
which keeps the energy finite. 

Moreover, the corresponding barrier is basi- 
cally impenetrable for the longitudinal W . In- 
deed, in the semiclassical approximation one 
can see that the action / p{r)dr diverges. There 
are basically two sets of levels; inside and out- 
side of the bag. 



2. 1 + 1 dimensions 

The content of this section is not strictly 
speaking relevant to the discussion of gauge 
quanta binding to a bag. It serves only the pur- 
pose of illustrating what exactly happens with 
the longitudinal mode in the simplest case of 
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a very large bag (l>{x) possessing a zero (ap- 
proximated by a 1-dimensional kink-like Higgs 
profile). We will then return to the finite-size 
bag with a Higgs zero. 

The ID kink separates two vacua, with op- 
posite signs of the Higgs VEV. By symmetry, 
4>{x) is assumed to be an odd function, being 
zero at the location of the kink, which we take 
at X = 0. Therefore, the effective mass of SM 
particles is zero near a; = 0, therefore creating 
attraction to this point. Let us discuss gauge 
bosons moving only in x direction normal to the 
kink. Transverse polarization modes E, M are 
decoupled from the longitudinal one and satisfy 
the same equation of motion, namely. 



d^F{x) 
dx"^ 



+ M^(t){xfF{x) = uj^F{x){k20) 



The effective potential is thus M^(/)(r)^, just 
as it would be for a scalar field. If ~ x 
at small a;, this potential is quadratic, with 
oscillatory-like levels with positive squared fre- 
quency w^. 

The longitudinal polarization, however, con- 
tains several additional terms 



d'^L(x) 



dx"^ 



+ MlrC^{xfL{x) 



2 dLdcj) 
)(a;) dx dx 



2 d 



(j) dx'' (f)'^ dx 

The first derivative of the function may be re- 
moved by the following substitution 



L{x) 



(A21) 



leading to the Schrodinger-like equation for the 
new wave function ipix) 



uP'4'{x) 



d^tPix) 
dx^ 



+ M^^ixf) (A22) 



dx^ 4>{x) 



ct^ixf 



Note that the new terms in the effective po- 
tential are singular at zeros of (j){x). The term 
with the second derivative of the kink shape 
4>{x) is harmless because a generic expansion of 
an odd function 4){x) = Cix + Csx^ + ... with 
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FIG. 9: The effective potential for the W levels 
V/{WOGeV) versus the distance x{WOGeV) from 
the kink center. The lower (green) curves are for 
the M modes, the upper (red) ones are for the 
longitudinal L ones. Three cases shown are for 
A/{100GeV) = 0.1,0.45 and 1, from bottom to 
top, respectively. 



some finite constants Ci , C3 will only produce 
a constant term in the effective potential, equal 
to 6C3/C1. Yet the other term has a very sig- 
nificant singularity. 

Now, let us take for definiteness a particular 
kink profile, e.g.. 



(x) = tanh{Ax) 



(A23) 



and see what modification the new term in- 
troduces to the effective potential. The effect 
depends crucially on the kink's width, or pa- 
rameter A. This is demonstrated in Fig. [9j 
as one can see in the lowest pair of curves, 
for A/{10QGeV) — 0.1, case there remains 
some attractive part of the kink, while in the 
A/(100GeV^) — 1 Upper) case there remains 
only repulsive potential, with the middle case 
of A/{100GeV) = 0.45 being the critical one. 
Apparently, the longitudinal mode is well above 
the transverse modes, if it even has any bound 
states at all. 

The reader can also be reminded about a 
comparison with the well known situation with 
fermions in a kink field, see recent discussion 
and original refs in [TU]. In this case the ef- 
fective potential also develops a negative part. 
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induced by the spin terms. However in this case 
the lowest level has exactly zero value. These 
correspond to the well-known fermionic zero 
modes, enforced by topological index theorems. 

Summarizing the lessons from the one dimen- 
sional kink problem: The lowest fermionic lev- 
els are in this case at zero energy, lower than 
the zero-point energy for transverse W-boson 
levels. The longitudinal W-bosons are strongly 
repelled from the zeros of the Higgs and thus 
their energy levels are much higher than those 
of transverse modes. 



b. Hamiltonian 
The Dirac Hamiltonian is given by 

Hd= [ d^r^P^hoi^, 



with 



ho =7° {—il''di + m. 



(B6) 



(B7) 



Therefore, 



Appendix B: Fermions in a Higgs 
Background 

1. 3+1 dimensions 



a. Spherical coordinates 

We take signature (H ) and use the stan- 
dard representation for Dirac matrices 

and ^ = ip'^^^ . In two component notation 



In spherical coordinates, we write 

(B3) 

with normalization J dr [F^ + G"^) — 1 and 
where ^jim are spherical 2-component spinors 
satisfying 



(B4) 



and 



doflij,n^l'j'm' ^ Sii'6jj'S„i,n>. (B5) 



Hr 



d^r 



'iif{<T ■ d)x + h.c.^ 
+ mcl){Tj*rj-x*x) ■ (B8) 



Using the spherical spinor properties one can 
show that 

{a-p)X = -^ (g'-^g), (B9) 
r \ r / 

(and similarly for rj) which leads to 



Hn^ dr 



-G' 



-G)f+(f' + -F]G 
r / \ r 



-m6 F^ - G^ 



(BIO) 



If ijja denotes a spinor solution with positive en- 
ergy Ea, i.e. ho tpa = Sai'a thcu, usiug ( [BQ] ), the 
equations of motion for F and G read 

{ea-mc^)Fa = -Gl+'^Ga, (Bll) 

r 

{ea + m4>)Ga = F^+'^Fa. (B12) 
r 



It is readily seen from (BIO) that using Eqs 



(Bll B12|, the total fermionic energy is again 



given simply by 



(B13) 



where ria is the occupation number of each state 
and N = ria is the total fermion number. 
Due to charge conjugation symmetry, we shall 
consider Uq > and > for both quarks and 
anti-quarks. 
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2. 1+1 dimensions 

We take signature (H — ) and use the standard 
representation 



7 



2 2 1 

CT , 7 = cr . 



(B14) 



for Dirac matrices in 1+1 and ip = V'^f^. We 
write a Dirac spinor as 




with left and right moving spinors 



(B15) 



1 



(F + G) 




The Hamiltonian reads 

Hd= J dx [F'G - FG' + m<j){F'^ - G^)] . 

(B16) 

Note that the 3d system becomes one dimen- 
sional in the r — ^ cx) limit, with the indices 
(j7m) becoming "internal" and F and G be- 
coming the upper and lower components of a 
1+1 Dirac spinor. 



3. Numerical results 

In order to see what actually happens with 
the fermionic levels as the bag's node moves 
from large distances towards the origin, we have 
numerically solved the Dirac equation in the 
background of the Higgs ansatz (42). 



In Fig{TO]we follow the I — Dirac wave func- 
tions for a kink profile. When the node is lo- 
cated at a rather large distance from the origin, 
the solution is very close to the analytic Gaus- 
sian solution discussed above: the maxima of 
both F{r) and G(r) are very close to the po- 
sition of the Higgs' node and \F\ ~ |G|. As 
the Higgs' node moves closer to the origin, al- 
though F{r) and G(r) are still approximately 
Gaussian, their maxima no longer trace the po- 
sition of the zero and they are visibly differ- 
ent, \F\ > |G|. This tendency continues as the 





-0.5 



FIG. 10: Upper and lower components (upper 
and lower solid lines) of Dirac's equation in the 
Higgs background (dashed line). The upper fig- 
ure shows a large kink profile with the lowest level, 
K = —1, asymptotically becoming a left-moving 
(i.e. F = — G) zero mode. As the node moves 
towards the origin, the system becomes increas- 
ingly non-relativistic and |-F| > |G|, as expected. 
A corresponding behavior for an anti-kink profile 
also occurs, now with the k = 1 level asymptoti- 
cally becoming a right-moving (i.e. F — G) zero 
mode. 



node moves past the origin and finally disap- 
pears. This is also easy to understand: as the 
bag potential becomes more shallow the system 
becomes increasingly non-relativistic, in which 
case one would indeed expect |F| :s> \G\. Apart 
from the I — level, it is also interesting to fol- 
low the behavior of the energy levels with I > 
as a function of R. 
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